MAGNETIC CONTRIBUTION TO THE SPECIFIC HEAT OF Pbi_^Eu^Te 
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The temperature dependence of the magnetic specific heat has been studied experimentally and 
theoretically in the semimagnetic semiconductor Pbi-xEuajTe for a;=0.027 and a;=0.073, over the 
temperature range from 0.5 K to 10 K, in magnetic fields up to 2 T. There was a maximum in 
the magnetic specific heat between 1 and 3 K even in zero and low magnetic fields; this maximum 
shifted toward higher temperatures with increasing magnetic field. The experimental data have 
been analyzed in the framework of a model in which we assume that the ground states of europium 
ions are split even without an external magnetic field. We present arguments which support this 
assumption and we show that it is possible to find a physical mechanism leading to the splitting 
which can explain the experimental results. 

PACS numbers: 71.70.Ch,75.40.Cx 



I. INTRODUCTION 

Semimagnetic semiconductors (SMSs), also known as 
diluted magnetic semiconductors (DMS), have been stud- 
ied extensively during the past two decades. Recently, 
there has been a considerable interest in these materi- 
als because of their possible applications in spintronics. 
Optical and magnetic measurements have shown that, 
in general, IV- VI SMSs with a 3d element as the mag- 
netic ion have a much weaker exchange interaction than 
that found in II- VI SMSs with the same magnetic ions. 
Also, the IV- VI SMSs with rare earths have a weaker 
exchange interaction than that found in the same mate- 
rials with a 3d element as the magnetic ion (for a review 
see Refs. Electron spin resonance investiga- 

tions in SnTe with Mn, Eu, and Gd have shown that the 
exchange interaction between free carriers and magnetic 
ions is roughly an order of magnitude smaller for Eu and 
Gd than for Mui^ The mechanism of the exchange inter- 
action among magnetic ions in IV- VI SMSs is still not 
well understood. 

Our previous investigations of the magnetic proper- 
ties of Pbi_2;Mn2;Te and Pbi_a:Eua:Te indicated a small 
{J/ks < 1 K), antiferromagnetic exchange interaction 
among magnetic ionsifi*2^ In Pbi_a;Eua;Te the absolute 
value of the exchange constant was about three times 
smaller than in Pbi_j;Mnj;Te and decreased with the 
increasing Eu-content. In Ref. ter Haar et al. have 
observed magnetization steps in the high-field magne- 
tization at mihkelvin temperatures in Pbi_a;Eua;Te and 
found exchange constant values similar to ours. By com- 
parison with the results in II- VI SMSs we came to a 
conclusion, that in IV- VI SMSs the dominant exchange 
mechanism is the superexchange between nearest neigh- 
bors (NN). In order to develop a more complete model 
and to obtain parameters for the exchange interaction, 
we have made complementary measurements of the mag- 
netic specific heat of Pbi_a;Mna;Te and analyzed the 



results together with the results of the magnetization 
and magnetic susceptibility measurements pifi It turned 
out that the mechanism of the exchange interaction in 
Pbi_2;Mna;Te may be more complex than just the NN su- 
perexchange. To explain the temperature and magnetic 
field dependence of the specific heat of Pbi_a;Mna;Te it 
was necessary to take into account a splitting of the 
ground-energy state of single Mn-ions in Pbi-ajMn^jTe 
and the p — d coupling between magnetic ion spins and 
free carriers. 

In the present paper we report studies of the mag- 
netic specific heat of Pbi_a;Eua;Te crystals and compare 
the results with those obtained in Pbi_a;Mna;Te. Some 
preliminary data have been recently reported!^ This is 
a first investigation of the magnetic contribution to the 
specific heat in rare-earth-doped SMSs. In the following 
sections we present the experimental results and analysis 
of specific heat measurements and give a physical model 
for the magnetic specific heat of Pbi_a;Eua;Te. 

In theoretical analysis presented in Section III we ar- 
gue that the experimentally observed magnetic specific 
heat is mainly due to single europium ions split in the 
disordered crystal environment. The splittings of mag- 
netic ions caused by the crystal field are small and have 
been usually detected in EPR experiments. However, for 
europium in IV- VI semimagnetic semiconductors these 
splittings are large enough to be observed also in dif- 
ferent kinds of measurements. In closely related semi- 
conductors Pbi_a;Eua;S and Pbi_a:Eua:Se magnetization 
steps due to the splitting of single europium ions have 
been observedii^iii In Pbi_a;Eua;Te we expect a similar 
effect. 

There is no consensus in the literature concerning 
the mechanism of the ground state splitting of rare 
earth ions in crystals. Several models were proposed and 
analyzed li^ii^ii^ It seems that at present it is impossible 
to decide unequivocally which mechanism should be ap- 
plied to a specific case. Instead, we suspect that in each 
case several mechanisms should be considered. There- 
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fore, in Section III we consider four different physical 
mechanisms leading to the splitting and estimate magni- 
tudes of the resulting splittings. 

The important part of the theoretical analysis is the in- 
corporation of deformations of the Pbi_a;Eua;Te crystal. 
These deformations are caused by the difference between 
Eu and Pb atoms and it turns out that even relatively 
small departures from octahedral symmetry may lead to 
splittings of the order 1 - 10 K. Therefore, they should be 
taken into account in analysis of magnetic specific heat 
measurements. 

The theoretical considerations in the present paper are 
limited mainly to the case of the a;=0.027 sample. For 
this sample most of the magnetic ions are singles, i.e. 
they have no nearest magnetic neighbors. For the sam- 
ple with higher concentration of Eu, a;=0.073, although 
the theory may be applied formally, its quantitative pre- 
dictions are not very reliable, because for such a high 
concentration a significant number of the Eu atoms is in 
larger, many atom clusters and the theoretical analysis 
of such a system is much more difficult. 



II. EXPERIMENT 

We have measured the specific heat of Pbi_a;Eua;Te 
with X values of 0.027 and 0.073. The samples of 
Pbi-j^Euj^Te were grown by the Bridgman technique and 
the Eu concentration was estimated from the amounts of 
the components introduced into the growth chamber and 
measured by energy dispersive x-ray analysis (ED AX). 
The nominal x-values were 0.03 and 0.06 with uncer- 
tainty of about 20%. The crystals were cut in the shape of 
Hall bars with typical dimensions 1.5 x 2 x 6 mm^. The 
samples were p type with carrier concentrations, from 
Hall measurements, of about 1 x 10^* cm~'^. With in- 
creasing X the hole concentration decreased and the mo- 
bility increased. 

Previously we have measured high-temperature mag- 
netic susceptibility and low-temperature, high-field mag- 
netization of Pbi_a;Eu2;Te with x up to O.li^S By fit- 
ting the susceptibility data to the Curie- Weiss law we 
obtained the average Eu content in our samples (xav) and 
very small Curie- Weiss temperatures indicating an anti- 
ferromagnetic exchange between Eu ions, J/ks = -0.38 K 
and -0.27 K for Xav = 0.027 and 0.073, respectively. 

The measurements of the heat capacity were performed 
in a cryostat using "^He and ^He systems over the temper- 
ature range 0.5 - 15 K in magnetic fields 0, 0.5, 1, and 2 
T. We used the standard adiabatic heat-pulse methodjii 
Errors in the heat capacity values were about 5%. The 
experimental details have been described elsewherei^li 

In order to obtain the magnetic contribution to the 
specific heat, Ch, it was necessary to subtract the spe- 
cific heat of the Pbi_a;Eua;Te lattice from the measured 
total specific heat of Pbi_j;Eu2;Te. This was not a sim- 
ple task. Bevolo et al. found that the specific heat of 
PbTe has an anomaly below 5 K and could not be fit- 



ted with the standard expression C = 7T -I- aT^, where 
7T and aT^ are the electronic and lattice contributions, 
respectivelyii^ In fact, they could not obtain a satisfac- 
tory fit to their data with an expression of the form, 
C = -fT + aT^ + X;r=i 5iT^'^^ unless n was at least 10. 
Therefore, we measured the heat capacity of our own 
Bridgman-grown PbTe sample in zero magnetic field and 
2 T, over the temperature range from 0.5 to 15 K and 
found that the temperature dependence was the same for 
and 2 T within our experimental error (as expected). 
In Ref. we described this experiment and have shown 
the specific heat of PbTe vs temperature in zero magnetic 
field. In our preliminary paper we have shown the result 
of simple subtraction of the PbTe specific heat from the 
Pbi_a;Eua;Te specific heatpii In the present paper we take 
into account the effect that the replacement of Pb with 
an atomic mass of 207.2 by Eu with an atomic mass of 
151.97 leads to a decrease in heat capacity, even for small 
values of X. To account for this we divided the entire set 
of PbTe specific heat data by empirically determined fac- 
tors, 1.11 for a;=0.073 and 1.04 for a:=0.027, before sub- 
tracting from the Pbi_j;Eu2;Te. These factors were de- 
termined by assuming that at temperatures above 15 K, 
in the absence of an applied magnetic field, the magnetic 
contribution to the specific heat of Pbi_£cEu£cTe is negli- 
gible. Therefore, this division by 1.11 (1.04) gave results 
for PbTe that were the same as those for Pbi_a;Eua;Te 
at 15 K for x = 0.073 (0.027). Since this is an empirical 
correction, we emphasize in the present work the data at 
temperatures below 5 K where the lattice specific heat is 
much smaller than the total specific heat. In the inter- 
esting region, below 2 K, the specific heat of PbTe was 
more than 3 orders of magnitude smaller than that of 
Pbi_:rEu:rTe. 

The magnetic specific heat data for Pbi_a;Eua:Te are 
shown in Figs.HandlS We believe that the scatter in the 
data represents the experimental error. For both values 
there is a broad maximum in the magnetic specific heat at 
about 2 K in zero magnetic field. The maximum is several 
times higher than that predicted by the cluster model 
of superexchange interaction between nearest neighbors. 
At higher magnetic fields the value of specific heat at the 
maximum increases and above 0.5 T it shifts to higher 
temperatures; for x = 0.073 the shift is smaller than for 
X = 0.027. This behavior is different from that observed 
in Pbi-^Muj^Teia or Sni-^Mn^^jTeyiS, where the value of 
the magnetic specific heat at zero magnetic field and all 
higher fields was nearly the same, but the shift of the 
maximum with increasing magnetic field was bigger than 
in Pbi_^Eu^Te. 



III. THEORETICAL ANALYSIS 

Before we start the analysis of magnetic specific heat 
(MSH) let us recall the model of the Europium atom 
in Pbi_a;Eua;Te which describes quite well magnetization 
and magnetic susceptibility experimental datai^ 
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the Hamiltonian for the spin subsystem reads 
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FIG. 1: Magnetic specific heat of Pbi_:i,Eu:rTe with x = 0.073 
in various magnetic fields. 
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FIG. 2: Magnetic specific heat of Pbi-:cEua:Te with x = 0.027 
in various magnetic fields. Points - experimental data, lines 
- theoretical predictions of the nearest neighbours interaction 
model for B=0 (continuous line) and B=0.5 T (broken line). 



The electron configuration of a free Eu atom is 
4/^5s^5p^6s^. It is believed that when it replaces a Pb 
atom in PbTe the electrons from the outermost shell, 6s^, 
play the role of electrons of Pb and contribute to crys- 
tal bindings. Due to the strong relativistic downward 
shift, the energy position of 6s^ Pb electrons is deeply 
in the valence bandiSfi As a result we obtain Eu^"*" ion, 
electrically inactive with respect to the crystal, despite 
the fact that it replaces Pb atom from IV group. Such 
a picture is confirmed by the fact that the presence of 
even 10% of Eu atoms has almost no effect on carrier 
concentrationjSi 

According to Hund's rule the ground state of Eu ion is 
^S, it means that seven 4/ electrons form the spin S=7/2 
and the total angular momentum L=0. Thus, the ground 
state of the ion is eightfold degenerate. This degeneracy 
is removed by external magnetic field, B or by interac- 
tion with another magnetic ion. With such assumptions 
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where the g-factor g = 2, /i^ is the Bohr magneton, and 
Jij is the exchange integral between the i-th and j-th 
spins. If the content of Eu is small, of the order of 1- 
3%, we may safely assume that most of the ions have 
no nearest magnetic neighbors and only small percent- 
age of them form two or three-atom clusters called, in 
the literature, pairs, open triangles, and closed triangles. 
Assuming the statistical distribution of Eu atoms in the 
PbTe lattice, one knows the average number of singles 
and the average number of atoms in pairs, open trian- 
gles, and closed trianglesjSi Then every thermodynamic 
quantity, in particular the magnetic specific heat, may be 
calculated. 

Although such a cluster model is successful in descrip- 
tion of magnetization in Pbi_a;Eua;Te^ it fails in the case 
of MSH. In Fig. |2 we see that the calculated MSH is 
much smaller than the one observed experimentally. The 
calculated MSH in Fig. |2l for B=0 is due to pairs and 
triples only, because the contribution from singles, which 
are for B — eight fold degenerate, is zero. The spe- 
cific heat due to larger clusters has, for a sample with 
a;=0.027, a very small contribution of about 2.5 %. We 
also think that theories like extended nearest neighbor 
pair approximation (ENNPA)22i based on the long range 
mechanism of spin - spin interaction, are not applicable 
to Pbi_a;Eua:Te bccausc even the nearest neighbor Eu-Eu 
exchange integral in Pbi_a;Eua:Te is small, J/ks «-0.25 
K^iS, and more distant interactions, which quickly decay 
with the distance, cannot explain the broad maximum for 
MSH for B=0. Let us also note, that for nonzero mag- 
netic field the theoretical curve above 1 K lies well below 
the experimental points. This picture gives evidence for 
a non-neghgible density of energy states in the energy re- 
gion far above 2 K, which must be taken into account to 
describe the experiment properly. This density of states 
of the system does not result from the model described 
by the Hamiltonian Eq. (QJ. 

The situation is somewhat similar to the case of 
Pbi_a;Mna;Te, where the peak in zero magnetic field 
MSH is also unexpectedly highiiS There is, however, at 
least one important difference between Pbi_a;Eua;Te and 
Pbi_a;Mna;Te. According to statistical mechanics, for 
any physical model describing MSH of a spin system, 
in particular for a model described by Eq. (QJ, we have 
following entropy relation: 



dT 



Ch{T) 
T 



^xR\n{2S+l)-kBHgH), (2) 



where R is the molar gas constant, S is the magnetic ion 
spin (for Eu^"*" S = 7/2), Ch{T) is the molar magnetic 
specific heat at temperature T and in an external mag- 
netic field, and qh is the degeneracy of the ground state 
of the system. The degeneracy depends on the magnetic 
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FIG. 3: Magnetic specific heat in Pbi-^Mnr^^Te with x=0.056 
in two magnetic fields (Ref. \l(J^ . 



field. In an external magnetic field B 0, gH=i and 
the second term on the right hand side disappears; thus 
calculating the difference between the left hand side of 
Eq. IpJ for magnetic fields B ^ and B = we obtain 
information about the degeneracy of the ground state 
in zero external magnetic field. Estimations, based on 
experimental results presented in Ref. 0, suggest that 
in the case of Pbi_a;Mna;Te this difference is nearly zero 
(see Fig. El- Thus, in the case of Pbi_a;Mna;Te, any model 
leading to gH=o 1 must be rejected. This is not the 
case of Pbi_a;Eua;Te. In Fig.[21we see that the difference 
between MSH measured in B = and in i? = 0.5 T is 
much bigger in Fbi-rEu.rTe than in Pbi_a;Mna;Te. As- 
suming, like in Ref. [ifl the linear temperature depen- 
dence of MSH below 0.5 K we find that the value of 
the integral for B = 0.5 T is equal to 0.41 J/mole K. 
For a;=0.027 the right hand side of Eq. ^ returns the 
value 0.47 J/mole K. The agreement is quite good. For 
B = 0, with the same assumption concerning the be- 
havior of MSH at the lowest temperatures, the value of 
the integral is 0.2 J/mole K. Because the first term on 
the right hand side, xRhi(2S + 1), does not change, it 
means that the second one, fcBln((7//), for B — must 
be positive. Thus, the experimental data suggest that 
in the case of Pbi_a;Eua;Te in magnetic field B=0 we 
have non-negligible degeneracy of the ground state of 
the spin system. This is the important difference be- 
tween Pbi_a;Mna;Te and Pbi_a;Eua;Te. In Pbi-^^Mn^j^Te 
even in B=0 the ground state is non-degenerate. The 
lack of degeneracy in Pbi_a;Mna;Te has been discussed 
in Ref. Efl 

At this point some clarifying remarks are necessary. 
First, we do not claim that Pbi_a;Eua;Te is a system con- 
tradicting the third law of thermodynamics. In reality, if 
we take into account all interactions and system degrees 
of freedom, the ground state is non-degenerate. How- 
ever, in our description we limit considerations to the 
spin subsystem and Eq. Q is derived and may be ap- 



plied only to such a subsystem. Secondly, we have no 
data at temperatures below 0.5 K, a region which may 
significantly contribute to the value of the integral. That 
is why the estimations of the left hand side of Eq. 12 are 
only semiquantitative and cannot serve as a rigorous jus- 
tification of the approach introduced below. However, we 
think that they provide important insight into the prob- 
lem and to some degree confirm the considerations given 
below. 

In the present paper we propose that the experimen- 
tally observed MSH in zero magnetic field is due to the 
splitting of the energy levels of the single Eu ions. From 
the discussion of different aspects of splitting that are 
presented later, it turns out that this splitting is caused 
primarily by two mechanisms: the disordered crystal field 
potential, which leads to virtual 4/^ 4/^5(i^ transi- 
tions, and the internal spin - orbit coupling on 4/ shell 
in the excited, 4/^5(i^, state. 

According to the Kramers theorem, the ground state of 
a single ion (of a system consisting of seven, i.e. an odd 
number of electrons on 4/ shell) is at least two-fold degen- 
erate in the absence of an external magnetic field. Thus, 
the ground state of the spin subsystem, a significant part 
of which consists of such split noninteracting ions, is also 
degenerate. Such an approach is in accordance with the 
experimentally observed difference between the integrals, 
Eq. (01, for zero and nonzero magnetic fields. On the 
other hand, the split singles contribute to the magnetic 
specific heat; therefore we expect our model describing 
the experiment to be better than the calculated curves 
in Fig. 12 In the following analysis an important role is 
played by 5d levels of europium. Let us discuss now the 
origin and the meaning of these states involved in the 
excited configuration 4/^5d^ of an ion. 

In the investigations of Pbi_2;Eu2;Te by Krenn et al^ 
the optical dipole transitions from the 4/ level of Eu 
to the vicinity of the bottom of conduction band were 
studied. The photon energy of these transitions for Eu 
concentrations corresponding to ours was found to be less 
than 1 eV. Notice that this kind of transition may take 
place only between states of different parity. Because the 
/ functions have odd parity, the states near the bottom 
of the conduction band must contain states of even 
parity. In pure PbTe the wave functions of the bottom of 
the conduction band, of Lg symmetry, are of odd parity. 
Thus, to explain the existence of the optical transition 
we must assume that in Pbi_a;Eua;Te the wave functions 
in the vicinity of the bottom of conduction band have 
certain components of even parity. The question arises 
about the origin of these components. The most natural 
candidates are 5d levels of Eu. Since it is known that 
in EuTe the conduction band is built mainly from 5d 
and 6s states of europium^^i one may expect that the 
addition of a few percent of Eu to PbTe will result in a 
contribution of 5d states to the conduction band states. 
Of course, the presence of even parity states may be 
also related to the disorder introduced by addition of 
Eu atoms to PbTe because, strictly speaking, the group 
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theory considerations apply only to the perfect crystals 
and for Pbi_a;Eua;Te containing several percent of Eu 
their conclusions cannot be taken too rigorously. In the 
approach presented below we model these states by a 
single, localized level of d symmetry, which we refer to 
as the 5d level of Eu. 



A. Crystal field potential in disordered crystal 

In our approach, the virtual 4/^ 4/^5c?^ transitions 
are caused by the electrostatic crystal-field potential. In 
the present subsection we estimate the order of magni- 
tude for this quantity in a disordered Pbi_2;Eua;Te crys- 
tal. 

In the crystal field theory a single 4/ or 5d electron 
moves in a potential that may be expanded into the fol- 
lowing infinite series: 



1=0 m=-l 



(3) 



where ro «0.5 A is the atomic length unit, 
and the functions Cim{d,(p) are related to spheri- 
cal harmonics Yim{0,yy) by the relation Cim.iO,(p) — 

/ \ 1/2 

( 2rfT ) ^;m(^', v)- In our approach we use the sim- 
plest model of crystal field potential, it means we assume 
that the crystal field potential is due to six point charges. 



each of which has charge Ze, placed at ViOiipi 



1,...,6 



with respect to the europium ion. Then the coefficients 
Aim readSi 



A, 



Ze 



6 
1=1 



(4) 



where C;*^ means complex conjugate to Cim- In our cal- 
culations we assume that Z — 2. 

Due to the difference between Eu and Pb atoms, also 
the Eu-Te and Pb-Te distances in Pbi-ajEu^jTe are dif- 
ferent. This difference causes local lattice deformations. 
These deformations are not limited to the nearest neigh- 
borhood, but extend over larger distances (several lat- 
tice constants) . If the concentration of Eu atoms is very 
small, the average distance between Eu atoms is large. 
Then, although the crystal lattice is locally deformed, 
the symmetry of the Eu surrounding is still preserved. 
However, with the increasing europium content, the de- 
formations originating from the different atoms start to 
overlap. Due to a random placement of Eu atoms in 
the lattice, we expect a random deviation of Eu-Te bond 
orientations from those in the perfect crystal. It turns 
out that these deviations cause significant ground state 
splittings of Eu^"*" ions. 

Let us consider a Pbi_a;Eu;rTe crystal with Eu content 
of several percent. To estimate the order of magnitude 
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FIG. 4: (a) Distribution of bonds' deviations from the ideal 
crystallographic directions for two different Eu compositions, 
(b) Distribution of Eu-Te distances in units of lattice constant 
of PbTe ao = 6.46 A. 



of the bonds' defiections we performed a numerical sim- 
ulation. First, we modeled a perfect lattice of crystalline 
PbTe containing 50'^ unit cells of PbTe. Next, a certain 
percentage of randomly chosen Pb atoms were replaced 
by Eu atoms. The cations were connected to anions by 
springs with equilibrium distances cJeu-Tc = 3.3 A and 
dpb-To = 3.23 A. Due to the lack of experimental data, 
the bond length dEu-To in Pbi_2;Eua;Te was taken as half 
of EuTe lattice constant. After applying zero tempera- 
ture Monte Carlo procedure, the equilibrium configura- 
tion of the lattice and the deviations of Eu-Te bonds from 
the perfect crystallographic directions for every Eu atom 
have been found. The typical bond deviation is of the 
order of several degrees. In Fig. E^a) we plot the prob- 
ability distribution for these deviations for two different 
europium contents. As it may be expected, the average 
deviation increases with the Eu content x. In Fig. ^b) 
we plot the distribution for Eu-Te distances. As we see 
this distribution is very well localized around the average 
distance. The relative changes of the bond lengths are of 
the order of 0.1%. 

The above, purely mechanical, model of disorder serves 
only to estimate the order of magnitude of deviation due 
to the difference between Eu-Te and Pb-Te bond lengths. 
It neglects, for example, the difference in strength of the 
bonds or the angular forces. In addition IV- VI com- 
pounds are very often disordered due to the presence of 
cation vacancies or granular structure of the material. 
Due to these reasons we expect bond deviations to be 
larger than estimated from Fig. 31 

In the description of our experimental data the degree 
of disorder will serve as one of the fitting parameters. 
We describe now the model of disorder which we apply 
to this fitting procedure. 

Using a Gaussian random number generator, we gen- 
erate random deviations of bond directions for each of 
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FIG. 5: Comparison of probability distributions for bond 
deviations (a) and crystal field potential (b). The continuous 
lines are the results of Monte Carlo simulations for a sample 
containing a;=0.027 of Eu and the broken lines are calculated 
using model of disorder described in the text for (^o=2.5°. 



the Eu atoms. More precisely, for a given Eu-Te bond, 
for example the one along (100) direction in the perfect 
crystal, we generate two random angles (} and 7, both 
with zero mean and the standard deviation equal to 4>{). 
It means that the considered bond will be characterized 
in spherical coordinate system by angles = 7r/2 + /3 
and (j) = 1- With slight modifications, a similar proce- 
dure is applied to the remaining five bonds of the given 
Eu atom. In our model of disorder we neglect changes 
of Eu-Te distances, because the simulations suggest that 
these changes are very small. For the configuration of Te 
atoms obtained in this way we calculate the crystal field 
potential, Eqs. 

In Fig. we compare the results of Monte Carlo proce- 
dure for a sample a;=0.027 with those resulting from the 
model of disorder introduced above for (^0 = 2.5°. We see 
that for such value of 4>q the differences between proba- 
bility distributions for one of the crystal field coefficients, 
^10, are not drastic. 

From the above analysis we see that the presence of dif- 
ferent cations leads to a deformation of the Pbi_a;Eua;Te 
lattice. In particular, deformations of nearest neighbour- 
hoods of Eu ions cause splittings of Eu ions ground states. 
In the next subsection we discuss physical mechanisms 
leading to the splittings. 



B. Mechanisms of ground state splitting of Eu 

ion 

In this subsection we estimate the magnitudes of 
the ground state splitting due to different physical 
mechanisms 



1. Direct influence of the crystal field on 4f electrons. 

Due to the strong spin orbit coupling for 4/ electrons, 
the Eu^"*" ion is not in a pure state, but the higher 
energy states of the 4/*" configuration are admixediA as 
shown: 



\S) = afSy2) 



7/2> 



cfD 



7/2/ 



(5) 



Thus we see that the total angular momentum L ^ Q 
and the ion may interact with the crystal field. The val- 
ues of the coefficients a, b, c for europium may be read 
from Table 5.1 of Ref. ^ a=0.986, 6=0.167, and c=- 
0.011. For such a state, using tables of spectroscopic 
coefHcientSfSS it is possible to derive matrix elements for 
{SM\Vcr\SM')^ For example, a term in Eq. © pro- 
portional to A20 yields the following term, H'^j^,j, in the 
effective spin Hamiltonian for 4/^ shelli^*2i 



2^/5 
105 



21 

T 



bcA2o{{r/rQf), (6) 



where M is the projection of the total spin on a quan- 
tization axis, M — —7/2, ...7/2, and the average of 
the second power of radius for 4/ shell wave function 
for europium, in atomic units, is ((r/ro)^)=0.938 (see 
Ref. 113) • 

Let us define A as the energy difference be- 
tween the highest and the lowest energy levels of the 
split ion. Then for the Hamiltonian Eq. A — 

(24y(5)/105)&c((r/ro)^)A2o. From simulations we know 
that the average value of A20 is of the order of 0.01 eV. 
Thus A is of the order of 10~^ meV which corresponds to 
0.1 K. Other terms in the expansion, Eq. (pj give split- 
tings of the same order of magnitude or smaller. In order 
to explain the magnetic specific heat we need A/fc^ to 
be of the order of 1 - 10 K. Thus the mechanism based 
on the direct infiuence of the crystal field potential on 4/ 
electrons cannot explain the splitting responsible for the 
experimentally observed magnetic specific heat. 



2. 4/^ ^ band states hybridization 

In 1978 Barnes et al^ noticed that the hybridization 
between the 4/ shell and the band states leads to a split- 
ting of the ground state of rare earth ions. The main 
idea of the model is to consider the excited states of the 
system in which the number of electrons on an ion's 4/ 
shell changes by ±1. According to the Hund's rule, in 
the excited states, 4/^ and 4/^, the total angular mo- 
mentum of 4/ electrons is nonzero. Taking into account 
internal spin - orbit coupling, the authors of Ref. ^3 ob- 
tained an effective spin - lattice interaction leading to the 
ground state splitting. Let us concentrate in this Section 
on processes which lead to 4/^ ^ 4/^ transitions. These 
processes may be important for Pbi_2;Euj;Te because, ac- 
cording to the optical measurements performed by Krenn 
et al^, the energy ei necessary to transfer an electron 
from the 4/ shell to the conduction band is of the order 
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of 0.5 eV, which is rather small. The model has been 
described in Ref. and has also been re-derived for the 
octahedral symmetry of the ion's neighborhood as dis- 
cussed in Ref. For completeness we describe it very 
briefly emphasizing the differences necessary to account 
for disorder in the crystal. 

The ground state of the system is a Eu ion in 4/^, 
conflguration plus the Fermi sea of electrons. This eight- 
fold degenerate state of the system is characterized by 
-7/2 < M < 7/2, where M is the projection of spin 
7/2 on a quantization axis which we take along the (001) 
crystallographic direction. In the excited states we have 
the ion in the 4/^ conflguration plus one additional elec- 
tron above the Fermi level, which is characterized by the 
set of quantum numbers q. This set of quantum numbers 
contains the wave vector from the flrst Brillouin zone, the 
number of the band, and the additional quantum number 
necessary to fully characterize the band state. This addi- 
tional quantum number enumerates Kramers conjugated 
states. (For semiconductors, for which the band spin- 
orbit coupling may be neglected, i.e. spin of band carrier 
is a good quantum number, one may think of this addi- 
tional quantum number as the projection of electron spin 
on a quantization axis. This is not, however, the case for 
Fbi-^EusTe where the band spin-orbit coupling cannot 
be neglected.) If we assume the validity of Hund's rule 
for the 4/® conflguration, L=3 and S=3, the Hamiltonian 
for the ion in the excited state reads: 



4/6 



A4yL • S. 



The hybridization elements are: 

{L,S,q\H\ M) 



(7) 



(8) 



The state \LzSzq) is the excited state of the system in 
which the projection on the quantization axis of the total 
angular momentum and spin of the ion are Lz and Sz, 
respectively, and there is one additional electron char- 
acterized by q above the Fermi energy. The coefficient 
(^—\^l^+i^ t/'2+2(j M^ may be derived using explicit forms 
of antisymmetric many electron functions for ion's states 
\LzSz) and | M). The element {q\h\(j)-L^a) describes hy- 
bridization between band state q and the 4/ spin orbital 
4>-L,a- The band wave functions are calculated within 
the tight binding model2i and the hybridization elements 
between 4/ shell and Te 6p and 6s orbitals are described 
by three constants V^/^, Vpf^ and T4/<t- According to 
Refs. Island 113 



Vpfcr — Vpfcr 



{rpr'f) 



1/2 



TOO 



d5 



^/tT — VpfT 



1/2 



Too 



d5 



(9) 
(10) 



where rjpfa- — 10\/2T/7r, rjpfj^ = — ISa/T/^/tt, rp=15.9 A, 
r/ =0.413 A and d is Eu-Te distance. We assume that 



^s/cr = "^p/o" the results of the calculations do not 
depend crucially on this assumption. The dependence of 
interatomic matrix elements on the direction of the Eu-Te 
bond with respect to the crystallographic axes of the ideal 
crystal is calculated according the method proposed in 
Ref. 13 The other details of calculations of the effective 
spin Hamiltonian are presented in Ref. 0. 

Due to the strong localization of the 4/ shell, the inter- 
atomic matrix elements responsible for the 4/ shell band 
states hybridization are very small, of the order of 0.1 eV. 
This is the main reason that there is a small splitting A 
despite the smallness of the transfer energy ei which is 
approximately 0.5 eV. The average splitting A/ks is of 
the order of 0.1 K. This is again much too small a value 
to explain the magnetic specific heat. 



3. 5d 



band states hybridization 



Unlike the 4/ orbitals of the Eu ion, the 5d orbitals 
are much more extended in space. Thus we expect the 
overlapping with neighboring Te orbitals to be much big- 
ger. Recently one of the authors (A. L.) has proposed 
a new mechanism leading to the ground state splitting 
of rare earth ions in which the 5d shell of a rare earth 
ion provides a bridge between the 4/ electrons and the 
rest of the crystalfi^ More precisely, an electron from the 
valence band jumps virtually onto the 5d level of the rare 
earth ion and interacts via a Heisenberg type of exchange 
with the spin of the 4/ shell. The Hamiltonian for eu- 
ropium in this excited state, 4/^5(i^ is of the following 
form 



- J/rfS • s + AsrfL • s, 



(11) 



where the first term describes the exchange interaction 
between 4/ spin S and the spin s of the 5d electron. The 
second term describes the spin orbit interaction on the 
5d shell. In the more general treatment of the problem in 
Ref. we used two spin orbit constants. Here we use a 
single spin orbit constant, X^d, and neglect the influence 
of the crystal fleld on 5d electron. For the disordered 
neighborhoods of europium atom which we consider in 
the present paper these simpliflcations do not change sig- 
niflcantly the flnal results. Like the previous subsection 
the influence of the surrounding comes via the hybridiza- 
tion between the 5d and the band states. This hybridiza- 
tion is described by three interatomic Eu-Te matrix el- 
ements Vpda, Vpd-K, and Vsda for which we assume fol- 
lowing Eu-Te distance dependenceSi Vpda ~ ^°dcr('^o/c')^; 
Vpi. = VX{ao/2d)\ and Vsda = V^iao/dy/^. Here 
ao=6.46A is the lattice constant of PbTe and d is the ac- 
tual Eu-Te distance along the given bond. The assumed 
values of three constants T^"^^=-1.5 eV, ¥^^^=0.7 eV 
and Vp^^=-1.6 eV are close to the ones used in Ref. 
in calculations of EuTe band structurei^ The values 
J/d=0.2 eV and A5d=0.08 eV have been taken from Ta- 
ble III of Ref.H3- The important parameter of the theory 
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is €2, the energy necessary to transfer an electron from 
the top of the valence band to the 5d level. Contrary to 
the case of gadolinium in PbTe where €2 is of the order 
of 0.4 eV, for europium we only know that it should be 
largeitiS. In our calculations we have assumed £2=! eV. 

For the above values of parameters the average split- 
ting A/ks is larger than for the two previous mecha- 
nisms, of the order of 0.5 - 1 K. It remains too small, 
however, to explain the magnetic specific heat. 

4. 4f ^ 4f^5d^ transitions 

The last mechanism of the splitting, which we con- 
sider, is based on 4/^ 4f^5d^ virtual transitions. As 
it has been already discussed earlier we treat 5d states 
not as the pure atomic states but hybridized with the 
band states. In some sense this mechanism is comple- 
mentary to the one considered in subsection 2 because in 
calculations in subsection 2 the Pbi_a;Eua;Te band states 
have not contained 5d orbitals. In other words, we may 
say that we add a certain amount of 5d to the 4/ states. 

The Hamiltonian of the model is constructed as fol- 
lows. In the ground state of the ion there are seven 
electrons on the 4/ shell. The total angular momen- 
tum is zero, the total spin equals 7/2 and this state is 
eight fold degenerate. As in previous sections |M) de- 
notes a state of the ion where the projection of the total 
spin on a quantization axis in 4/'' configuration is M 
{M = -7/2,. ..7/2). 

The excited state configuration is 4/^5d^. It is described 
by following Hamiltonian: 

Hipbd^ = Hip + Asdl • s - J/dS • s -I- Vcr + eo, (12) 
where 

7?4/6 - A4/L-S + A:i^(L-S)' (13) 

is a Hamiltonian for six electrons on 4/ shell describing 
splitting of F state due to 4/ spin orbit interaction. We 
assume that in the excited configuration, six electrons 
on 4/ shell behave according to the Hund's rule, i.e. 
their total spin S'=3 and the total angular momentum 
L = 3. The values of the spin orbit coupling constants 
A4/=0.03 eV and A^y=0.0005 eV have been fitted to de- 
scribe properly a splitting of the 4/^ configuration cal- 
culated from first principles (see Table VIII in Ref. 1^3) . 
The second term in Eq. lfT2|l describes spin orbit cou- 
plings on the 5d shell. The next two terms correspond 
to exchange interactions between the 4/ and bd spins 
and the crystal field potential acting on electron on the 
bd shell. The last term, eo, is the energy necessary to 
perform the 4/^ Af^'ad^ transition. We neglect the 
infiuence of the crystal field potential on the 4/ electrons 
since we have checked that its infiuence on the final result 
is very small. 

The crystal field potential, Eq. (EJ, enters the Hamil- 
tonian of the model in two places: in the term describing 



the excited states, Eq. H12|l . and in the terms describing 
4/^ ^ 4/^5d^ transitions. The basis in which we de- 
scribe excited states of the ion is denoted by \SzLzlzc), 
where Sz and Lz correspond to projections on a quantiza- 
tion axis of total spin and total angular momentum of six 
electrons on 4/ shell, respectively, while lz and cr are z-ih 
components of angular momentum and spin of the sev- 
enth, the bd electron. Using properly antisymmetrized 
wave functions or the concept of fractional parentage 
coefHcients^SiSi, we have derived the following form for 
the hybridization elements: 

{SzLzlz<J\Y.l=lycr{r^)\M) ^ (14) 

where and (j)^/ are bd and 4/ orbitals, respectively. 
The external magnetic field B is taken into account by 
adding to the Hamiltonian Zeeman term: 

HB=g^iBB{Sz + <Jz) + ^^BB{Lz + lz), (15) 

where g-factor g=2 and is the Bohr magneton. Calcu- 
lating the matrix elements of the Hamiltonian in the basis 
|M), \SzLzlzCr) we obtain a 498x498 matrix. The eigen- 
values of this matrix enable us to calculate the magnetic 
specific heat according to the standard rules of statisti- 
cal mechanics. The first eight eigenvalues of the Hamil- 
tonian matrix correspond to the split levels of the Eu^+ 
ion. The higher levels describe the 4/^5ci^ configuration 
and they are separated from the lowest eight ones by an 
energy of the order of eo. That is why the splitting A 
is defined here as the difference between eighth and the 
lowest eigenvalue. The value of the matrix element of r 
between radial 4/ and 5d wave functions, necessary to 
calculate hybridization elements, may be estimated from 
the Table VI of Ref. ]^ In calculations we assume that 
(4/|r/ro|5c?)=l. The contribution from terms in Eq. (pj 
with I > 1 may be omitted since the coefficients Aim 
decay quickly with I and most important contribution 
comes from terms with 1=1. The other parameters of 
the model are the same as those in the previous sections, 
J/d=0.2 eV, A5d=0.08 eV, and eo=l eV. For such values 
of parameters of the model we obtain A/fc^ of the order 
of 1-5 K. 

From the analysis of different mechanisms of the split- 
ting we conclude that the last one, i.e. the one based 
on 4/^ ^ 4/^5(i^ transitions leads to the largest ground 
state splitting. This mechanism, with the above values of 
parameters will be used to explain the magnetic specific 
heat for a sample containing x=0.027 of europium. 

C. Magnetic specific heat 

The magnetic specific heat is calculated according to 
standard rules of statistical mechanics. We take into 
account singles, pairs and triples. We assume that the 
Hamiltonian for pairs and triples is of the form given in 
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FIG. 6: Magnetic specific heat for sample x=0.027. Points 
- experimental results. Continuous lines - theoretical calcula- 
tions for the mechanism described in Section III.B.4. 



Eq. Q with Eu-Eu exchange integral J//cb=-0.25 KiM 
For these clusters we neglect the splitting caused by the 
crystal field. We do not expect that the approximation 
would introduce a large error, because a total splitting of 
the pairs' energy levels in zero magnetic field is 28 • 2 J, 
in our case 28-2-0. 25^14 K; this is the same order of 
magnitude as the splitting of singles caused by the dis- 
ordered crystal field. The magnetic specific heat due to 
singles is calculated in the following way: First we gen- 
erate 100 random Eu environments as was described in 
Section III A. For each set of tellurium positions we cal- 
culate the energy levels of the split ion and we calculate 
the corresponding magnetic specific heat. As the mag- 
netic specific heat due to singles we take the average over 
these 100 samples. We found that 100 samples are suf- 
ficient; for more samples we obtained nearly the same 
result. The parameter (po is treated as the fitting param- 
eter. The best results have been obtained for 0o=3° and 
they are shown in Fig. El 



D. Discussion 

In the theoretical analysis of MSH we have concen- 
trated on the sample containing x=0.027 of Eu atoms. 
In this sample, if one assumes a random distribution of 
ions, more than 97% of them are in single, pair, and 
triple clusters and about 72% are singles. As we see from 
Fig. El the theoretical curves describe the experimental 
data quite well. The differences are probably due to the 
fact that the assumed model of disorder does not fully 
refiect all the complexity in a real crystal. This theoret- 
ical description has been achieved using only two fitting 
parameters: eo and (po. All other parameters are known 
from the literature. 

In the sample containing a;=0.073 Eu, only about 40% 



of ions are singles and more than 24% are in clusters con- 
taining more than three atoms. That is why the quan- 
titative analysis is more difficult. However, the experi- 
mental results presented in Fig. are in semiquantitative 
accordance with the proposed model. 

Applying Eq. Q for B ^ and for B = we see that 
the difference of the left hand sides for these two cases 
equals 



dT 



T 



dT 



Ch=o{T) 
T 



ks ln(gH=o) 



(16) 

where gH=o denotes the degeneracy of the ground state 
of the system in zero magnetic field. (As in the previous 
analysis the ground state of the spin system in the pres- 
ence of magnetic field is nondegenerate, i.e. In^gn^o) = 
0.) Calculating the integrals using experimental data 
with the assumptions discussed previously (temperature 
dependence of the specific heat for r<0.5 K is linear) 
we obtain this difference equal to 0.17 J/mole K. For 
a sample containing 0.073 of europium 40% of Eu ions 
are singles. It means that one mole of Pbi_a;Eua;Te con- 
tains 0.4 • X ■ Na of europium singles, where Na is the 
Avogadro number. Other europium ions are in larger 
clusters. According to our model the ground state of 
each Eu single is doubly degenerate. Thus the degener- 
acy of the ground state of the spin system due to singles 
is 2°-^'^'^-** and the right hand side of Eq. Ijlfill equals 
OA - X ■ i?ln2. For x = 0.073 we obtain 0.16 J/mole K. 
Assuming that the ground state of Eu ions in larger clus- 
ters is non-degenerate we obtain a very good agreement 
with experimental data. However, we realize that the 
integrals in Eq. ltT6ll are estimated from a limited set of 
data. Therefore this agreement confirms our model only 
semiquantitatively. 

The important problem we should consider is whether 
the proposed model is consistent with the earlier magne- 
tization and magnetic susceptibility measurements. In 
the earlier theoretical description of the experimental 
magnetization data, using a Brillouin-function analysis 
for singles, we assumed that the Eu ions in zero mag- 
netic field are not split 

Fig.Hshows the magnetic-field dependence of the aver- 
age spin calculated with the proposed model, continuous 
lines, for two different temperatures. The broken lines 
show the behavior of Brillouin functions calculated for 
spin 8=7/2. We see that although at the lowest tempera- 
tures the differences are noticeable, they are small. Thus 
our analysis has shown that the magnetic specific heat 
measurements reveal properties of the system (density 
of states) which are not refiected in the magnetization 
measurements. 



IV. CONCLUSIONS 

We have measured magnetic specific heat of 
Pbi_2;Eu2;Te for a;=0.027 and a;=0.073. We have shown 
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that the experimental results may be explained assum- 
ing that the single Eu ions are split in a disordered crys- 
tal field potential even without external magnetic field. 
Analyzing the possible mechanisms of the splitting we 
have concluded that the main contribution to the split- 
ting comes from the virtual 4/"^ <-^- Af^bd^ transitions. 
We have also shown that the model is suitable for de- 
scription of the earlier measurements of magnetization. 
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FIG. 7: Magnetic field dependence of average spin of a single 
ion calculated according to the model - continuous lines. The 
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